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Abstract 

We compute the momentum-transfer dependence of the proton Panli form factor F2 in 
the endpoint overlap model. We find the model correctly reproduces the scaling of the ratio 
of F2 with the Dirac Form factor Fi observed at the Jefferson Laboratory. The calculation 
uses the leading-power, leading twist Dirac structure of the quark light-cone wave function, 
and the same endpoint dependence previously determined from the Dirac form factor Fi. 
There are no parameters and no adjustable functions in the endpoint model’s prediction 
for F2. The model’s predicted ratio F2(Q^)/Fi{Cp) is quite insensitive to the endpoint wave 
function, which explains why the observed ratio scales like 1 /Q down to rather low momentum 
transfers. The endpoint model appears to be the only comprehensive model consistent with 
all form factor information as well as reproducing fixed-angle proton-proton scattering at 
large momentum transfer. Any one of the processes is capable of predicting the others. 


1 Introduction 

The electromagnetic form factors know as Fi and F2 are an important probe of the internal 
structure of nucleons. A popular theoretical model assumes that at high momentum transfer 
these quantities can be factorized into a hard scattering contribution and a so-called distribution 
amplitude. The distribution amplitude has no information about the proton wave function except 
the parton momentum fraction Feynman-a: dependence and some spin factors of a short-distance 
expansion. The focus of the short-distance (SD) model [IJ|^ is a perturbatively calculable hard 
scattering kernel. The model generates an order by order expansion in powers of the inverse 
momentum transfer-squared, l/Q^. The expansion has often been claimed to be the unique 
prediction of QCD. However the task of comparing the model to the larger theory of QCD was 
never completed, and obviously cannot be explored within the SD model itself. 
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Yet model predictions can be compared to experimental data. The SD model predicts that 
F 2 {Q^)/Fi{Q^) ^ 1 /Q 2 for large Q 00' A simple way to obtain this uses a perturbative quark 
mass to flip two quark helicities in the internal lines. The experimental results obtained at the 
Jefferson lab [^[^, however, showed that the ratio F 2 {Q^)/Fi{Q'^) ~ 1/Q in the energy range 
2 GeV^ < < 5.6 GeV^. This contradicted the prediction of the SD model that some thought 

had been established. The results played an important role in dramatizing the failure of the SD 
model, which had also been anticipated earlier 10 It is now clear that the SD model might 
apply only at very large energies, which are inaccessible experimentally. Even at asymptotic 
energies there is no proof the model dominates. 

Since the SD model fails it is imperative to explore alternatives. 

Work by Miller et al, Lin et al, and Gloet et al 12 -14 has reproduced the experimentally 
observed momentum dependence of F 2 . These calculations emphasize the importance of the 
quark wave functions, i.e. the role of hadron structure, as opposed to the role of perturbation 
theory. 

Kivel and Vanderhaeghen 15 used Soft Gollinear Effective Theory (SCET) to analyze soft 
spectator contributions which represent a class of diagrams which also give the SD scaling of 
Q^T' 2 /Fi ~ const, for large Q^. It is possible that at smaller of order few GeV^, the soft 
spectator contributions might lead to the observed experimental behavior. However, it is not clear 
how to extract these contributions systematically. In an earlier analysis, Belitsky et.al 16 obtained 
the dependence F 2 {Q^)/Fi{Q'^) —>• l/Q^ * logiQ "^/which matches with the observed data, 
by introducing higher twist light cone wave functions. The logarithmic term is a result of an 
integration over the soft endpoint region, where the assumptions of the SD model no longer hold. 
Hence we find that some studies in the past 15 16 have attributed the observed experimental 


behavior of F 2 /F 1 to the contributions arising from the soft spectator quarks in the end point 
region. 

The significance of the end point region for the calculation of the ratio F 2 /F 1 is rather in¬ 
teresting in view of the recent claim that an end point model {EP model) can comprehensively 
explain the scaling behavior of many exclusive processes [T 20 ^ . The model relates the observed 
scaling to the behavior of the quark wave function as Eeynman-a: —>-1. In this limit one of the 
quarks carries most of the proton longitudinal momentum. The model appeared several times in 
the literature, yet it was dismissed prematurely, often for reasons that its premises contradicted 
the assumptions of the SD model. Eor that reason the EP region was long regarded as a nuisance. 
Many efforts attempted to show the EP model’s contribution would be suppressed, but the efforts 
were unsuccessful. 

Once given fair consideration, the EP model appears to provide the simplest explanation of 
several experimental observations. In [17] , we applied the EP model to compute the pion form 
factor, the proton Dirac form factor Fi((5^), and the proton-proton elastic scattering cross section 
at high momentum transfer. We found that one consistent wave function for the end point region 
could be extracted by fitting the experimental form factor data. The same wave function then 
predicts the scaling behavior observed in proton-proton fixed-angle scattering. We extend this 
study here in order to determine the proton Pauli form factor F 2 . We find that the formalism 
predicts F 2 {Q^) without introducing any new parameters. 

Let us briefly explain the physics. It is well known that quark mass insertions produce a quark 
helicity flip, which can ultimately produce the proton helicity (more specifically, chirality) flip 
characterizing E 2 . Quark mass terms are negligible in the high energy limit of the SD model. 
This is because they compete with terms scaling like the large momentum Q. The role of a quark 
mass is qualitatively different in the EP model. The soft quarks with momentum fractions x ^ 0 
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already have very small momenta. Their momenta are of the order of the QCD chiral symmetry 
breaking scale A. Then an equally small contribution from a quark mass is not a relatively small 
effect, and it cannot be neglected. Let us repeat that the attempt to banish small momentum 
regions from QCD never worked out. An unexpected consequence of small momenta appearing at 
leading power order is that mass effects can appear at the same order. 

Another effect makes this even more interesting. Under a Lorentz transformation with rapidity 
y in the z direction, the big light cone + component transforms by and the small component 


like 


All previous calculations known to us at leading power order integrate quark wave 


functions over the small momentum in the first step. This appears to be much more safe than 
integrating over the transverse momentum components, which scale like 1 compared to e“^. Yet 
we have discovered a limit-interchange error occurs. Integrating away the small components is the 
first step of the SD model producing a visible factorization into separated hadronic parts. The 
assumption, actually a hope, that some factorization dominates is what demands that step. Yet 
that step instantly causes F 2 to scale no larger than l/Q®. When the small momenta components 
are retained in the scattering process we find a contribution to F 2 scaling like 1/Q^. The integrals 
cannot be represented by effective, pre-integrated quantities that depend only on Feynman-x. 
This phenomenon contradicts the tenets of factorization. In the EP model, the leading power 
contribution to F 2 comes from an inseparable union of initial and final state proton states. 

Finally all of this occurs with one simple wave function, which happens to be the most often 
cited, leading twist example. There is no particular reason to favor leading twist coming from a 
short distance expansion. There is every reason to use a wave function of leading power in the 
large momentum P. It is seldom noticed that the leading power, leading twist wave function has 
both chirally-even and chirally-odd components. A single wave function can both maintain the 
proton’s chirality in Fi, and flip the chirality in F 2 . 

In section we show that by respecting the necessary integration region, while using the 
endpoint dependence of the proton wave function obtained in 17 , we obtain the experimentally 


observed scaling behavior for F 2 /F 1 . This is a remarkable prediction of the model: If attention 
had been given 30 years ago, it would have predicted F 2 in advance of the data. Reversing 
the argument, the observed scaling dependence of F 2 /F 1 predicts Fi and pp scattering at high 
momentum transfer. None of these facts requires appealing to an unusually large logarithmic 
correction, or an unusually large dimensionful scale. As far as we know it is the first time that 
one model is actually consistent with the known data. 

Quark orbital angular momentum is a topic of great interest. No orbital angular momentum 
(OAM) enters the SD model, because a theoretical preference for factorization demands integrat¬ 
ing over quark transverse momenta before the actual reaction has even been set up. Information 
about transverse size is lost by that step. When OAM is re-cast into a twist expansion 16 the 
sequence of operations dictated by the SD model produces a 1/Q® dependence for F 2 . Refer¬ 
ences [18[|19| showed that avoiding the SD assumptions and performing the transverse momentum 
integrations to compute F 2 led to power law dependence for F 2 intermediate between 1/Q'^ and 
1/Q®. That is, the integration region assumed to dominate asymptotically was not the actually 
dominant region, whether or not an endpoint issue was considered. While the asymmetry of the 
endpoint integration regions produces a rather obvious role for OAM, that is not the focus of 
this paper. This paper is about using the same leading twist Dirac and endpoint structure found 
in the Fi calculation to calculate F 2 . The calculation is relatively simple, and agrees remarkably 
with data. Even more remarkably, the ratio F 2 {Q‘^)/Fi{Q'^) is quite insensitive to the endpoint 
wave function, explaining why the observed ratio goes like 1/Q down to rather small momentum 
transfer. 
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2 Endpoint Calculation 





Figure 1: The basic kinematics of the end point contribution to the proton form factor. The photon with 
momentum q scatters with one of the quarks which carries the dominant fraction of the proton 
momentum. 


Here we describe the calculation of through the quark mass contribution. One quark is 
struck by the virtual photon. The remaining quarks will be in a small momentum region, such 
that their incoming and outgoing momenta are entirely determined by their wave functions. No 
interactions are computed for those particles, because perturbative interactions would double¬ 
count what is already included in the wave functions. We will use the same wave functions to 
compute F 2 as previously determined 


17 from Fi, found to be consistent with pp scattering. 


2.1 Coordinates 

The basic diagram for proton electromagnetic form factor is given by Fig. The initial and 
final proton 4-momenta are P and P, with q = P — P . Initial quark momenta kj (masses rrij) 
are unprimed, while final momenta use the same label with a prime. We let ki denote the struck 
quark, and ^ 2 , fca denote the spectators. Our coordinates are {energy, p^, Py, Pz)- We use a 
Lorentz frame where the incoming and outgoing protons momenta are 




p'r- _ 




2 Q Q Q 


g'^ = (0,Q,0,0). 


( 1 ) 


Here mp is the mass of the proton. 
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We introduce a basis for transverse momenta: 


=(0, 0, 1, 0) = y'; P ■ y = P' ■ y' = 0; 
n'^=^(0,-1, 0, -1), P-n = 0; 

n"^=^(0, 1, 0,-1), P'-n' = 0. 


Here P = P' = vectors along the direction of 

propagation of the incoming and outgoing protons respectively. The components of the quark 
3-momenta are expressed as 


h 

K 


Xi ^ P T kiyy — ( XiQ j‘2^^ 0, XiQ -t- ( kin / 2, kiy , kin /5 

v2 

x'^P +k[nn' + Kyy = (x'Q/2, 0, a;-Q/2)-h (fc'„/v^, k'iy, 


( 2 ) 


The four momenta of the quarks are then given by, 


kt = 



f j 0 Q kin 1 Q kin 
(k^,-Xi--^,hy,X,--^ 



^in y y ^ 

72 ’ *2 72 


( 3 ) 


2.2 The Matrix element 

With the electromagnetic current operator and N standing for Dirac spinors, the matrix 
element for the interaction is parameterized by 


< p's'l J^|ps >= —ie 


Pi(Q7(ivV^) 


^{Q^)N'ia^''qnN 


( 4 ) 


Let be the Bethe-Salpeter 3-quark wave function in the proton with spinor indices shown. 
Let symbol stand for the quark-photon vertex, propagator factors, and momentum conserva¬ 
tion factors, displayed in a moment. The model for the reaction is 


< p s'\j^\ps >= / n 


(Pki (Pk' 


\5^{ki + *2 + fca - P)5^{k'i + 7 + 7 - P') 


(27r)4 (274 
4'a'/3'7'(7) ^ aji-fiki) 


( 5 ) 


Here is 


■Ad'" = - ie-f'^,j'^{ki + q- 7)(7 - m 2 )y{k 2 - - m 3 )y^S*{k 3 - 7) 

- ie7^,^7(fci +q- 7)(7 - m 2 )y^S‘^{k 2 - 7)(fe “ m 3 )yj‘^{k 3 - k' 3 ) 

- +q- 7)(7 - m 2 )yj^{k 2 - 7)(fe “ m 3 )y{k 3 - k' 3 ). (6) 

Note the delta functions (54(7 —7)A^(7~7) "''^hich explicitly enforce momentum conservation 
of spectator quarks. 
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The initial light cone coordinates are defined as 


k+ = k° + 


XiQ 


k~ = ky - 


Vf' 


Final state symbols have a prime. In literature, it is standard to use the co-ordinates k~ = k~p'^; 
K'^ = k'^/p'^\ K~=k~p~^]K~^ = k~^/p^ which is just parameterizing the light cone co-ordinates 
with the momenta , p . 


2.3 Integration 

It is generally assumed that wave functions 'I''(fc'), 'P(fci) of 4-momenta are peaked near the on-shell 
region. In that region, the actual wave function can be replaced by its integral over the small mo¬ 
mentum component k~, k~, producing the usual light cone wave function Y'{x[, ^li)) Y{xi, k_\_i) 

221 . 


< p'sV'^lps >= J n + I^t + 1^3 - + k±2 + kj_3) (7) 

5(ki+ -f + K + - l)d^{k'j_i + fc^2 + j dnjd{K'{- + k'^ + 4" - 

j 

M^{ki,ki){jY[dKY6{K^ +K 2 + K 3 -ml)'i/ai3j{ki)) 


The rest of the calculation cannot use the same approximation, because the delta-functions vary 
rapidly: Hence the process is indivisibly linked together by the integrations. 

The above expression for scattering kernel AI^ has an important dependence on and Kg , 


which cannot be overlooked. This is the point where our calculation begins to differ from previous 
ones. 

The basic problem is that for the soft spectator quarks it is not reasonable to assume that 
their four momentum square, fc^, is approximately zero. We expect k^ to be of the order of A^. 
In a constituent quark model, these quarks are assumed to be approximately on-mass-shell with 
masses of the order of few hundred MeV for the up and down quarks. In general the behavior of 
the quark propagator is expected to be more complicated and one can model its form by solving a 
truncated Schwinger-Dyson equation 23 -25 . For our purpose, it is adequate and self-consistent to 
assume that is dominated by the on-shell region and the k~,k ~ dependence can be replaced 
by the on-shell expression. 


rrig -I- kj 


( 8 ) 


As explained above, we assume that the mass of the slow spectator quarks is of the order 
of a few hundred MeV. On the other hand, the struck quark is a perturbative object. Treating 
it consistently uses a mass of order of a few MeV. We ignore the tiny and power-suppressed 
helicity-flip contributions from the struck quark. 

Now, doing a change of variables gives the standard form with 
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< ps'\J^\ps >= 


J + *2 + iCa - l)(5^(fc_Li + k±2 + k±3)6{x[ + + Xg - 1) 

SHk'x, + + fcla) (^a'/3'y k'^.)M^Y^0^ix,, 

( 9 ) 


The delta functions of Eq.([^ and Eq. § lead to the following conditions, 


X 3 = 

= 1 — xi — 

^ 2 ; 

*'3 = 

l-x\- X 2 ] 

^In 

= —k 2 n — 


kly = 

- ^ 2 y ^3y7 

k' 

^In 

— —k' — 

V . 

^3n5 

U - 

Kly - 

- —U _ u . 

^ 2 y ^3y’ 

kiy 

= k'iy \ k2y 

II 


^In 


■ 7); 

h' 

^In 


^ 2 n 

= -p={-X 

72 

2 ); 

U — 

^2n — 



■,X2=X2+0 



The light cone wave function Y of leading twist and leading power of large P is 26 27 
In 




16^2 


■{{fC)o,p(-f^N)^V + {f-i^C)^pN^A + ^K.P"C)„;3(7''75iV)^r}. (10) 


Here V^A,P are scalar functions of the quark momenta, N is the proton spinor, N,. the number 
of colors, C the charge conjugation operator, a^i, = |[ 7 ^, 7 jy], and /at is a normalization. This 
wave function was previously used to compute Fi, and is now being applied to compute Ej. 

It may come as a surprise that the same chirality structure creating Fi can predict E 2 . Fig. 
shows a cartoon of the chirality flow. Each term in the collection has been classified as 

chirally even or chirally odd depending on whether it conserves helicity (even, anti-commutes with 
75 ) or flips helicity (odd, commutes with 75 ). Since momentum conservation is trivial it is not 
shown. The chirality flow of the V,A,F terms are shown at the top. A typical combination of 
diagrams flipping the final state proton chirality is shown at the bottom. This diagram needs 
one ( 1 ) internal flip of low momentum spectator quark chirality, which appears as the closed loop 
with a mass insertion indicated by “X.” The cartoon shows how the Dirac algebra works without 
needing to do the algebra. 

Returning to Eq. inserting the wave function Eq. |10[ and extracting the terms which lead 
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Figure 2 : Chirality flow in the calculation of F2, indicated by arrows. Chirally-even vertices conserve he- 
licity and chirally odd ones flip helicity. The standard leading twist wave function contains both 
types, shown across the top. A typical combination of diagrams flipping the final state proton 
chirality is shown at the bottom. This diagram needs one ( 1 ) internal flip of low momentum 
spectator quark chirality, which appears as the closed loop with a mass insertion indicated by 
“X.” 


to F2 yields 

_, j r I 

^ F2 = J dkiydk2ydxidx2-^ x { 

- m2)a'c.{h - 

+ - 'm 2 )a'aih “ m 3 )p’l3i{aapP^C)apT*V 

+ - m2)p-pi{ap,,,P.,C)^pT*T 

+ [^'757m.(I^3 + 

+ ...} (11) 

The 1/Q^ factor after the integration measure comes from the Q dependence of 5(k^ — k^) = 
5{{kt +x,Q/^j2)-{k + + x^Q/^/2)). Evaluating the first two terms in the above expression and 
isolating the F2 contribution gives 

_ I i f 1 

N -a>^‘'q^NF2 ~ J dkiydk2ydxidx2-^ 

The other terms are similar. 

2.4 The endpoint wave function and F 2 

The leading power wave functions of Ref. [17| were determined in the endpoint region: 

V,A,T cx (1 — Xi)xie~^i'F . 


2m p 


8mp[{P ■ k2)m3 + {P ■ ^3)7712] T*V 


( 12 ) 
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The exponential dependence on the transverse momentum is a generic form that restricts the 
range of xi G (1 — 1) and X 2 G (0, ^). 

The dot products are 


P t IQ^ ^ 2 

P ■h = ki\l — +m^- Xi— . 


In terms of the light cone variables, this gives 


Pi. m ^ 2 Q" 

P-ki =[kZ + x*^)y — +m^p-Xi — 


^ Q \ Q ^ Q2 


It follows that 


F 2 (X I dkiydk2ydxidx2-„2 


1 


X 8mp[(A(3)m3 + (A(5 )to 2] (1 - aii)e ^ (1 - x^e ^ 

1 


dxidx 2 ^Q(l — xi)(l — xi)e a 2 e 
1 


■Q^ 


= —Q— A = _ 

Q2^Q3Q Q5 


(13) 


2.5 The ratio of form factors 

In our estimate of the form factor F 2 we used the wave function given in Eq. |12[ whose x 
dependence was determined by fitting the Dirac form factor, Fi. However it is easy to see that the 
ratio F 2 /F 1 is independent of the precise form of the wave function within the end point model. 
Consider a rather arbitrary wave function 

V,Arcx/(xl)e-'=-/^^ (14) 


This leads to the Dirac form factor 17 


Fi oc y dkiydk 2 ydxidx 2 -^[ 8 Q'^m 2 m 3 ] /(xi)e f{x'i)e 


Similarly the form factor F 2 becomes, 


F 2 oc y (ifciydA; 2 ydxidx 2 y^ 8 mp[(A( 5 )TO 3 + (A( 5 )m 2 ]/(xi)e /( 7 )e ■ 

Taking the ratio gives 

F 2 ^8TOp[(A(5)m3 + (AQ)m2] l 
Fi ^[8(52 to2TO3] Q 


(15) 


(16) 


(17) 


Thus the ratio of form factors in the endpoint model is independent of the precise form of the 
wave function. 









S. Dagaonkar et al. - The Dirac Form Factor Predicts ... 


10 



The JLAB data shows QF^jF^ ^ constant starting from as low as 2 GeV^. At such 
low values F\ differs significantly from its high-scaling behavior, which is observed to set in for 


> 5 GeV^ 


In the low regime a more complicated wave function is needed to fit the 


data. However Eq. 17 follows quite generally since the dependence on the wave function cancels 
out while taking the ratio. 


3 Soft gluon exchange 


It can be verified that addition of low momentum gluons in the interaction will not change the 
scaling behavior of the Pauli Form factor F^. Consider the simple case of 2 gluon exchange 
illustrated in Fig.[^. 

The matrix element for this diagram is 


[dki][dk[] 




+mi) i i 


■j(^ +mi) 

[(—f5s7p) 2 _ )\a'c 

Pf2 ''^2 


Ph - P91 Pg 2 


-j5s7A]/3'/3 


where 


[dki] = dxidki±S{xi +X2 + x^- l)(5^(fci_L + fc2J. + ^3 _l) 


Evaluating the traces and extracting the co-efficient of N ia^''q,yN we find 


iv'- AT IP fui ^uin 4 8im2ipfi ■ P)iN ia^^’'q^N)T*V 

N i(j^''q^NF2 ~ / [dk,] [dfc']g^e —-y-- -P -G —■ 


{p\ - mf){p\ - ttiDpI^pI^ 


Keeping only leading power term for the limit Q 3> A, dropping transverse momentum integrals 
of order the hadronic scale and substituting V,T from EqfT^ gives 


F 2 


dxidx2dx'idx'2 




A2(-(l-a;i)Q2)A2A2 


X (1 - xi)(l - x'l) 
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Each integral dx over an interval of length A/Q contributes a power of 1/Q. The integration 
of 1 — a: over 1 — K/Q < x < 1 contributes a power of 1/Q- It follows that 


F 2 oc J dxidx'dx 2 dx 2 {i — x') oc 


1 


Q xQ X ...Q QS ■ 
Thus the gluon exchanges do not change the leading power behavior. 


4 Conclusions 


As mentioned in the Introduction, if the EP model had been given adequate attention 30 years 
ago, a fit to the known \/Q'^ dependence of Fi would have then predicted F 2 /F 1 ~ l/Q at large 
Q, just as eventually observed. The calculation was never done, despite the model’s visibility after 


initial development by Drell, Yan, Feynman, and others. 29 ■ 31 


Between then and now came a period attempting to dispense with hadron structure in form 
factors, and replacing protons with perturbation theory, which revealed very little about hadron 
structure. We find that one simple pattern of an endpoint wave function, previously determined 
in Ref. [17| and going like 1 — a;, explains many independent experiments. The endpoint region 
produces the original and earliest quark-counting model 30 . For each spectator integration dx 


restricted to a; ;< A/Q an integral goes like A/Q. For each hard struck quark with 1 — A/Q ^ a: < 1 
an integral goes like A./Q. Thus three quarks leads to Fi ^ l/Q^- The leading twist Dirac 
structure, which has no room for orbital angular momentum, still allows a reversal of the proton’s 
chirality characterizing F 2 , and F 2 ^/Q^- These are not asymptotic limits, but generic results 
of power-counting that apply in the region Q >> A, namely Q > GeV. 

The fact that QF 2 {Q^)/Fi{Q^) is nearly constant with Q down to rather low values of mo¬ 
mentum transfer is now understood. At small Q the details of the endpoint wave function enter 
the calculation, and replacing dx ^ A/Q is not accurate. It is possible to fit that dependence 
from data for Fj rather trivially. However the integrations for F 2 are so nearly like those for Fi 
that the details of the wave function cancel out in the ratio F 2 /F 1 . The rule that F 2 /F 1 ~ 1/Q 
for Q >>GeV naturally extends itself into the region of Q ^ few GeV. When future experiments 
probe higher momentum transfers we are confident that QF 2 {Q^)/Fi{Q^) will remain constant, 
regardless of what might occur with the numerator and denominator. 
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